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GENERAL DECAY IN SOME TIMOSHENKO-TYPE SYSTEMS 
WITH THERMOELASTICITY SECOND SOUND 

MOHAMED ALI AYADI, AHMED BCHATNIA, MAKRAM HAMOUDA, 

AND SALIM MESSAOUDI 


Abstract. In this article, we consider a vibrating nonlinear Timoshenko system with 
thermoelasticity with second sound. We discuss the well-posedness and the regularity 
of Timoshenko solution using the semi-group theory. Moreover, we etablish an explicit 
and general decay results for a wide class of relaxating functions which depend on a 
stability number p. 
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1. Introduction and setting of the problem 

Beams represent the most common structural component found in civil and mechani¬ 
cal structures. Because of their ubiquity they are extensively studied, from an analytical 
viewpoint, in mechanics of materials. A widely used mathematical model for describ¬ 
ing the transverse vibrations of beams is based on Timoshenko beam theory TBT (or 
thick beam theory) developed by Timoshenko in the 1920’s. The TBT accounts for 
both the effect of rotational inertia and shear deformation that occur within a beam as 
it vibrates. These factors are neglected when applied to Euler-Bernoulli beam theory 
EBT (or thin beam theory), which is appropriate for beams with small cross-sectional 
dimensions compared to the length. In fact, a fundamental assumption in EBT is that 
cross sections remain plane and normal to the deformed longitudinal axis throughout 
deformation, while in TBT cross sections remain plane but do not remain normal to the 
deformed longitudinal axis as the shear deformation is taken into account. The cross 
section rotation from the reference to the current configuration is denoted by p> in both 
models. In the EB model, this is the same as the rotation of the longitudinal axis. In 
the Timoshenko model, the difference is used as measure of mean shear distortion. 
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In 1921, Timoshenko [28] gave the following system of coupled hyperbolic equations 

pu u = {K(u x - ip)) x , in (0, L) x 1R+, 

Ipiptt = ( ELp x ) x + K(u t - <p), in (0, L) x H+, 

together with boundary conditions of the form 

EI<p x |^=o =°> («*-<P)l*2o =0. 

as a simple model describing the transverse vibrations of a beam. Here t denotes the 
time variable and x is the space variable along the beam of length L, in its equilibrium 
configuration, u is the transverse displacement of the beam and ip is the rotation angle of 
the filament of the beam. The coefficients p, I p , E, I and K are respectively the density 
(the mass per unit length), the polar moment of inertia of a cross section, Young’s 
modulus of elasticity, the moment of inertia of a cross section, and the shear modulus. 

System a, with the above given boundary conditions, is conservative and the 
natural energy of the beam, given by 

£(t) = \ f (p\u t \ 2 + I p \tpt\ 2 + EI\tp x \ 2 + K\u x - ip\ 2 ) dx, 

1 Jo 

remains constant in time. 

Vibration has long been known for its capacity of disturbance, discomfort, damage 
and destruction. Since a long time, many researchers have been investigating ways to 
control this phenomenon. However, with the development of control theory for partial 
differential equations over the last few decades, it is not surprising that the issue of sta¬ 
bility and controllability of Timoshenko-type systems has received a great attention of 
many mathematicians. One effective method for vibration control is passive damping. 
Damping is most beneficial when used to reduce the amplitude of dynamic instabilities, 
or resonances, in a structure. 

Damping is the conversion of mechanical energy of a structure into thermal energy. 
A structure subject to oscillatory deformation contains a combination of kinetic and 
potential energy. 

A damping effect may be caused by applying the beam to internal or boundary 
frictional mechanisms. Depending of the nature of the beam’s material, a damping 
effect may be rotating beam. For Viscoelastic materials with long memory, some beams 
are characterized by possessing both viscous and elastic behavior. As a result of this 
behavior, some of the energy stored in a viscoelastic system is recovered upon removal 
of the load, and the remainder is dissipated in the form of heat. 

Kim and Renardy [7] considered m together with two boundary controls of the 
form 

c)v 

K<p{L,t) — K—(L,t) 

Elfu.t) 

and used the multiplier techniques to establish an exponential decay result for the 
natural energy of (ED- They also provided numerical estimates to the eigenvalues of 
the operator associated with system (ED- An analogous result was also established by 
Feng et al. [3j, where the stabilization of vibrations in a Timoshenko system was studied. 
Raposo et al. m studied (ED with homogeneous Dirichlet boundary conditions and 


= a—(L,t) Vt> 0, 

= Vt> 0, 
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two linear frictional dampings. Precisely, they looked into the following system 


{ piu tt - K(u x -tp) + u t = 0, in (0, L) x IR + , 

P 2 Vtt ~ btp xx + K(u x - <p) + <p t = 0, in (0, L) x E + , 
u(0,t) = u(L,t ) = <^(0, t) = ip(L,t) =0, t > 0 

and proved that the energy associated with (jl.2ft decays exponentially. Soufyane and 
Wehbe [27] showed that it is possible to stabilize uniformly (11.11) by using a unique 
locally distributed feedback. They considered 


(1.3) 


putt = ( K{u x - <p)) x , in (0, L) x IR + , 

< IpVtt = ( EIip x ) x + IC{u x -<p) - bipt , in (0, L) x IR+, 
u(0,f) = u(L,t ) = <^(0,t) = <p(L,t) = 0, t > 0, 


where b is a positive and continuous function, which satisfies 


b(x ) > b o > 0, Vie [ao, ai] C [0, L], 


In fact, they proved that the uniform stability of (|1.3I) holds if and only if the wave speeds 
are equal ; otherwise only the asymptotic stability has been proved. Rivera 

and Racke m obtained a similar result in a work, where the damping function b = b(x) 
is allowed to change sign. They also in treated m a nonlinear Timoshenko-type system 
of the form 

j Pl^Ptt &1 {‘Pxt VOa: 0, 

\ P2^tt ~ X(^x)x + (?2{Vx, Ip) + dipt = 0 , 

in a one-dimensional bounded domain. The dissipation here is through frictional damp¬ 
ing which is only in the equation for the rotation angle. The authors gave an alternative 
proof for a sufficient and necessary condition for exponential stability in the linear case 
and then proved a polynomial stability in general. Moreover, they investigated the 
global existence of small smooth solutions and exponential stability in the nonlinear 
case. 

Shi and Feng [21] used the frequency multiplier method to investigate a nonuniform 
Timoshenko beam and showed that, under some locally distributed controls, the vi¬ 
bration of the beam decays exponentially. The nonuniform Timoshenko beam has also 
been studied by Ammar-Khodja et al. [2] and a similar result to that in [24] has been 
established. 

Ammar-Khodja et al. [T] considered a linear Timoshenko-type system with memory 
of the form 


(1.4) 


PiVtt ~ Kppx + ip) x = 0, 

P 2 'P>tt - bip xx + Jo g(t - s)ip xx (s)ds + K(tp x + ip) = 0, 
< ip(x, 0 ) = ip 0 (x), pt{x,0) = y>i(x), 
ip{x, 0) = ip 0 (x), ip t (x, 0) = ipi{x) 

. ^(0 ,t) = = Ip(0,t) = ip(l,t) = 0, 


in (0, L ) x IR + , and proved, using the multiplier techniques, that the system is uniformly 
stable if and only if the wave speeds are equal j and g decays uniformly. More 

precisely, they proved an exponential decay if g decays in an exponential rate and 
polynomially if g decays in a polynomial rate. They also required some extra technical 
conditions on both g' and g" to obtain their results. This result has been later improved 
by Messaoudi and Mustafa m and Guesmia and Messaoudi [5], where the technical 






4 


AYADI, BCHATNIA, HAMOUDA, AND MESSAOUDI 


conditions on g" have been removed and those on g have been weakened. Also, Guesmia 
and Messaoudi [6] considered the following system 


(1.5) 


PWtt ~ K(ip x + VOx = 0, 

P 2 tptt - ki/j xx + fog(t - T)(a{x)i/> x (r)) x dT + K(ip x + ip) + b(x)h(ip t ) = 0, 
< (p(x, 0 ) = ipo(x), <pt(x,0) = <pi(x), 

ip(x,0) = ip 0 (x), ip t (x, 0 ) = ipi(x), 

, ¥>(0 ,t) = tp(l,t) = ip(0,t) = ip{l,t) = 0, 


in (0,1) x IR + . They proved under similar conditions on the relaxation function g , which 
are similar to those in j2j, and by assuming that 


a(x) + b(x) > p > 0, Mx € (0,1), 


an exponential stability for g decaying exponentially and h linear, and polynomial sta¬ 
bility when g decays polynomially and h is nonlinear. 

Concerning stabilization via heat effect, Rivera and Racke m investigated the fol¬ 
lowing system 

{ PWtt ~ = 0, in (0, L) x IR + , 

P 2 ^tt ~ bi’ xx + K{ip x + ill) + 7 9 X = 0, in (0, L) x 1R + , 

PsOt- K0 xx + '-fil’ xt = 0, in , (0, L) x IR+, 

where are functions of (x,t) model the transverse displacement of the beam, 

the rotation angle of the filament, and the difference temperature respectively. Under 
appropriate conditions on a, pi, b, K, 7, they proved several exponential decay results for 
the linearized system and non exponential stability result for the case of different wave 
speeds. 

Concerning Timoshenko systems of thermoelasticity with second sound, Messaoudi 
et al. pj] studied 

PWtt ~ cr(tp x ,^)x + PPt = 0, in (0, L ) x IR + , 

P 2 ^tt - b'lpxx + k(cp x + V’) + PQ X = 0, in (0, L) x 1R+, 

P3&t + 7Qx + Si/itx = 0, in (0, L) x IR+, 

ToQt + q + k9 x = 0, in (0, L) x 1R + , 

where <p = ip(x, t ) is the displacement vector, ip = ip(x, t ) is the rotation angle of the 
filament, 9 = 9(x, t ) is the temperature difference, q = q(x, t ) is the heat flux vector, p\, 
P 21 P 3 , b, k, 7, 5, k, p, tq are positive constants. The nonlinear function a is assumed 
to be sufficiently smooth and satisfy 


0V* (0,0) = cr^(0,0) = k, 


and 

< ^V > xV>x ( 0 ) 0 ) <T i p x ^i > (0, 0 ) — 0 . 

Several exponential decay results for both linear and nonlinear cases have been estab¬ 
lished in the presence of the extra frictional damping ppt ■ 

Fernandez Sare and Racke (3j considered 


( 1 . 6 ) 


PPPtt ~ k(<p x + ip) x = 0, 

P2i’tt - bip xx + k(<p x + ip) + 59 x = 0 , 

P3&t + 7 q x + Siptx = 0, 

rq t + q + k9 x = 0, 


in (0, L) x 1R + , 
in (0, L) x 1R + , 
in (0, L) x 1R+, 
in (0, L) x IR + , 


and showed that, in the absence of the extra frictional damping (p = 0), the coupling 
via Cattaneo’s law causes loss of the exponential decay usually obtained in the case of 
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coupling via Fourier’s law m- This surprising property holds even for systems with 

history of the form 

(1.7) 

piptt ~ k(ip x + ip) x = 0, in (0, L) x 1R+, 

( pii'tt - btpx x + k{p x + VO + / 0 +o ° g{s)i}) xx {., t - s)ds + 86 x = 0, in (0, L) x 1 R+, 

P3&t + + Hhx = 0, in (0, L) x 1R+, 

T Qt + q + k0 x = 0, in (0, L) x 1R+, 


Precisely, it has been shown that both systems (11.61) and (11.71) are no longer exponentially 
stable even for equal-wave speeds . However, no other rate of decay has been 

discussed. 

Very recently, Santos et al. [22] considered (11.61) and introduced a new stability number 

Pi \ (P2 _ Pi\ _ PiS 2 Pi 
kp 3 ) \b k) kbp 3 ’ 



and used the semi-group method to obtain exponential decay result for p = 0 and a 
polynomial decay for p ^ 0 . 

The boundary feedback of memory type has also been used by Santos [21] . He 
considered a Timoshenko system and showed that the presence of two feedbacks of 
memory type at a portion of the boundary stabilizes the system uniformly. He also 
obtained the rate of decay of the energy, which is exactly the rate of decay of the 
relaxation functions. This last result has been improved and generalized by Messaoudi 
and Soufyane [9]. For more results concerning well-posedness and controllability of 
Timoshenko systems, we refer the reader to mm, m, m and [25] ;26]. 

In this paper we consider the following Timoshenko system: 

( 1 . 8 ) 




pirn - k((p x + ip) x = o, 

P2iptt - bipxx + k(<p x + VO + S6 X + a(t)h(ipt ) = 0, 
p-.iOt + q x + Hxt = o, 
rq t + /3q + 0 X = 0, 

<Px(0,t) = ¥V(M) = "0(0,7) = = q{0,t) = q(l,t) = 0, 

<p(x, 0) = <po(x), <Pt(x, 0) = <pi(x), 

= V’o(^), ^t{x, o) = ipi(x), 

0(x,0) = 0q{x), q(x, 0 ) = q 0 (x), 


in (0,1) x 1R + , 
in (0,1) x 1R + , 
in (0,1) x 1R + , 
in (0,1) x 1R + , 

V t > 0, 

V x € ( 0 , 1 ), 

V x € ( 0 , 1 ), 

V x € ( 0 , 1 ), 


where, pi, p 2 , p 3 , b, k, 5, (3 are positive constants, ip = <p{x,t) is the displacement 
vector, 'V = ip(x,t) is the rotation angle of the filament, 9 = 9(x,t ) is the temperature 
difference and q = q(x,t ) is the heat flux vector. Also, a and h are two functions to be 
fixed later. 


Using (| 1. 8 D i . (I1.8IU and the boundary conditions (II.81b . we have 


dt 2 


f 

Jo 


ip(x, t)dx = 0 and 


d f 1 

— / 9(x, t)dx = 0 . 
at J o 


Consequently, we obtain 
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If we set 


and 


ip(x,t) = <p(x,t) - 


j " ( pi(x)dx S jt + J ipo(x)dx), 


9(x,t) = 9(x,t) — f 9o(x)dx , 

Jo 

then (<p,ijj,9,q) satisfy also the system (11.81) . and we have 

/ <p(x,t)dx = 0 and / 9(x,t)dx = 0 . 

Jo Jo 

From now on, we use the new variables (<p,ip,9,q), but we denote them by (ip,ijj,9,q), 
for simplicity. 


The article is organized as follows. First, in Section^ we use the semi-group theory 
to prove the existence and uniqueness of solutions of system (11.81) . Next, in Section [3l 
we study the asymptotic behavior of the energy of solutions of system (j 1 . 8 D using the 
multiplier method. For that purpose, we assume some hypotheses on a and h. The op¬ 
timal exponential and polynomial decay rate estimates can be obtained in some special 
cases with explicit nonlinear terms. 


2. Well-posedness and regularity 

In this section, we discuss the well-posedness of the problem (II.8j) . using the semi¬ 
group theory. We consider the following hypotheses on a and h: 

(Hi) : a : M+ —>• M+ is differentiable and decreasing. 

(A 2 ) : h : M — > M is a locally Lipschitz function satisfying h{ 0) = 0. 

We introduce the Hilbert space: 

L *(0,1) = {v € L 2 (0,1) : [ u(s)ds = 0}, 

Jo 

H l{o, 1 ) = if^o, 1 ) n lI(o, 1 ), 

#*( 0 , 1 ) = {u € H 2 (0, 1) : 'Ua.(O) = v x (l) = 0}. 

The energy associated with the system (11.81) is defined by: 

1 r 1 

E{<p, V*, 9, q)(t) = - / (pup* + P2^t +wl + Kvx +^f +P 3 O 2 +rq 2 )dx. 

1 Jo 

Let 

H = H\{ 0,1) x Ll( 0,1) x Hq (0,1) x L 2 (0,1) x L 2 (0,1) x L 2 (0,1), 

be the Hilbert space endowed with the inner product defined, for 
U = (ui,U 2 ,u 3 ,U 4 ,u 5 ,u 6 y € H, V = (vi,v 2 ,v 3 ,V 4 ,v 5 ,v 6 ) t € H, by 

(t 7 ) V) H = Pi (U 2 ,V 2 ) L 2 (0j1) + P 2 (u 4 , V A ) L 2 (0j1) + k (u lx + U 3 , V lx + V 3 ) L 2 ( 0)1 ) 

+b (u 3x , V 3x ) L 2 ( 0 , 1 ) + P3 ( u 5, v 5)l 2 (0,1) + r ( u Oi v o)L 2 (0,1) ' 
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For <f> = (ip, u, 0, v, 9, qf and <f > 0 = (y? 0 , <fii, 0o, 01, # 0 , Qof, where u = ip t and v = 0 t , 
()1.8p is equivalent to the abstract first order Cauchy problem 

r d -$(t) + {A + B)$(t) = 0 , VtGM+, 

$( 0 ) = $o, 

where A : D{A ) C H —> H is the linear operator defined by 

/ -a \ 


( 2 . 1 ) 


dt 


( 2 . 2 ) 


and B : D{B) C H 


A<£ = 


\ 


—u 

~~PxX — 0 ;X 

pi ' pi 

—v 

~^0 xx + ^dPx + 0 ) + -j^0 x 

~k qx + ~h Vx 


k + ¥ 

H is the nonlinear operator defined by 

( 0 \ 


= 


0 

0 

a{t)h{v ) 
0 
0 


The domain of the operator A is given by D{A) = {d> G H ; AQ £ H} and endowed 
with the graph norm 

II^IId(a) = \\®\\h + II^IIld 

can be characterized by 

D(A) = (H 2 ( 0,1) n HUo, 1)) x Hl{ 0,1) x (H\ 0,1) n Hl( 0,1)) 

X-f^o(O) 1) x Hl(0, 1) x Hq(0, 1). 

The domain of the operator B is given by D(B) = {d> G H ; B& G H} = H. 

We first state and prove the following lemmas which will be useful to deduce the well- 
posedness result. 

Lemma 2.1. For d> G D{A), we have (Al<f>,<l>)# > 0. 

Proof. For any $ = {ip, u, 0, v, 9, q)* G D(A ), we have 


(A<h, §) H = k -{u x + v){ip x + i/j))dx+ {—kifxx — kif x )udx + b -u x 0 a 

Jo Jo Jo 

+ / (-bfJxx + k(<Px + 0) + S9 x )vdx + {q x + 5v x )0dx + {/3q + 9 x )qdx. 

Jo Jo Jo 

Using integration by parts and the boundary conditions in (|1.8|) . we obtain 

(At<h, &)h = (3 [ q 2 dx > 0. 

Jo 

This ends the proof of the lemma. 


dx 


Lemma 2.2. I + A is a surjective operator. 
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Proof. For any W = (w\,W 2 ,w 3 ,W 4 ,w 3 ,wq) € H, we prove that there exists V 
(v 1 ,v 2 ,v 3 ,V 4 ,v 5 ,v 6 ) <E D(A) satisfying 

(/ + A)V = W. 

That is, 


(2.3) 


— V2 + V\ = Wl, 

-kv lxx - kv 3x + pivi = pi(wi + w 2 ), 

-V4 + V 3 = W 3 , 

~bv 3xx + k(v i x + v 3 ) + 5v 5x + p 2 V 4 = P 2 W 4 , 
v 6x + 5v 4x + p 3 v 5 = p 3 w 5 , 

{P + r)v 6 + v 5x = twq . 

Then (12.31) i. (12.31) ^ and (12.3k yield 

(2.4) v 2 = v 1 -w 1 €Hl(0,l), 

(2.5) V 4 = v 3 - w 3 € i?o(0,1), 

V6 X = P 3 W 5 + <5^30; - Sv 3x - p 3 v 5 . 

By integration over (0,x) and using va(0) = 103 ( 0 ) = 1 ) 3 ( 0 ) = 0, we obtain 


( 2 . 6 ) 


I 


vq = p 3 / w 5 ds + Sw 3 - 5v 3 -p 3 v 5 ds. 


f 

J 0 


We substitute (|2.6I) into CEIDe and we get 

px px 

p 3 / w 5 ds + 5w 3 - Sv 3 -p 3 v 5 ds 
. Jo Jo 


= TW 6 . 


V5x + (P + r) 

Hence, we deduce that 

px px 

(2.7) — v 3x +(ft+ t)5v 3 + p 3 (j3 + t) / v 5 ds = (P + r)5w 3 + (P + T)p 3 w 5 ds-rw 6 

Jo Jo 

Again, we substitute (12.71) into (J2.3l) .i . we get 


-bv 3xx + kv i x + kv 3 + 5 

-(P + r)p 3 [ 
Jo 


(P + t)5v 3 + p 3 (P + t) f v 5 ds - (P + t)5uj 3 

Jo 


w^ds — twq 


+ P2V 3 = P2{w 3 + W4), 


and we infer that 

px 

- bv 3xx + kv lx + kv 3 + 5 2 (P + t)5v 3 + p 3 5(P + r) v 5 ds + p 2 v 3 = (P + t)5 2 w 3 

Jo 

pX 

(2.8) +{P + r)5p 3 / w 5 ds - 5rw e + p 2 (w 3 + W 4 ). 

Jo 

By using (12.71) . (|2.8I) and (12.31) 9. it can be shown that v\, v 3 and 05 satisfy 

(2.9) 

-kv lxx - kv 3x + pivi = hi <E Ll( 0 , 1 ), 

~bv 3xx + kvi x + kv 3 + (5 2 (P + t) + p 2 )v 3 + p 3 S(P + r) / v 5 ds = h 2 € T 2 (0,1), 


px ^ 

-PoVbx + Po{P + t)8v 3 + p 2 (P + t) v 5 = h 3 € L 2 (0,1), 

^ Jo 
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where 


hi = pi(wi + w 2 ), 

h 2 = {/3 + t)5 2 w 3 + (/3 + t) dp 3 f* w 5 ds - 5 tw 6 + p 2 {w 3 + xu 4 ), 
h 3 = P 3 (P + t)5w 3 + (/3 + r)pl f Q x w 5 ds - p 3 rw 6 . 

Let u = ( 111 , 113 , 1 x 5 ) and v = (v\, v 3 , 1 x 5 ), a simple multiplication of (12.91) 1 . (1 2.9 D o and 

rx 

(12.91) -^ . by xxi, 1 x 3 and / u 3 ds respectively, and integration over (0,1) yield 

Jo 


(2.10) • —k I vi xx uidx — k I v 3x u\dx + p\ / viuidx = I hiuidx , 

Jo Jo Jo Jo 

fi fi fi r 1 

, / 1 1, / ,1,.. 1 1, / „_1 lx 2 

Jo 


/ ui^xxidx - / ' 

Jo Jo 

—6 / v 3xx u 3 dx + k / uia;?X 3 dx + k / v 3 u 3 dx + (<5 2 (/3 + r) + P 2 ) / v 3 u 3 dx 

Jo Jo Jo Jo 

r 1 /*1 

+/ 03 < 5(/3 + r) / ( / v 5 ds)u 3 dx = / h 2 u 3 dx, 

Jo Jo Jo 

r 1 rx r 1 /*# 

P3 / V 5 x( u 5 ds)dx + p 3 (/3 + t)5 / u 3 ( / u 5 ds)dx + 

Jo Jo Jo Jo 

/*1 rx rx r 1 rx 

pI(/5 + t) / (/ v 5 ds)( u 5 ds)dx = / h 3 ( u 5 ds)dx. 

Jo Jo Jo Jo Jo 


Using integration by parts and the boundary conditions yield 


S' 


k I vi x uixdx + k v 3 ui x dx + pi viuidx = / hiuidx, 

Jo Jo 


f 


•b v 3x u 3x dx + k / vi x u 3 dx + k / v 3 u 3 dx + (<5 2 (/3 + t) + p 2 ) / v 3 u 3 dx 
Jo Jo Jo Jo 

r 1 /*£ /*1 

H-/c>3<5( / S + t) / ( / v 5 ds)u 3 dx = / h 2 u 3 dx, 

Jo Jo Jo 

r 1 /*1 rx 

•P 3 v 5x u 5 dx + p 3 (/3 + t)5 / u 3 ( / u§ds)dx + 

Jo Jo Jo 

/*1 /‘Ol /* 1 C r 1 rx 

P 3 ((J + t) / (/ v 5 ds)( u 5 ds)dx = / /i 3 ( / u 5 ds)dx. 

Jo Jo Jo Jo Jo 

The sum of the previous equations gives the following variational formulation 
(2.11) b(v, u ) = l(u), 

for all u = (ui,u 3 ,u 3 ) € H}(0 ,1) x JJq(0 ,1) x L 2 (0,1), where b is defined by 

b(v, u) = k (yix + v 3 )(ui x + u 3 )dx + pi viuidx + b v 3x u 3x dx 

Jo Jo Jo 

r 1 r 1 rx r 1 

+(<5 2 (/3 + t) + p 2 )/ v 3 u 3 dx + p 3 5(/3 + T) (/ v 5 ds)u 3 dx + p 3 v 5x u 5 dx 
Jo Jo Jo Jo 

r 1 /*# /*1 rx rx 

+P3{/3 + t)6 v 3 ( u 5 ds)dx + pl(/3 + t )/ (/ v 5 ds)( u 5 ds)dx, 

Jo Jo Jo Jo Jo 

by 

/*1 r 1 /*1 rx 

l{u) = / hiuidx + / h 2 u 3 dx + / /z 3 ( / u 3 ds)dx. 

Jo Jo Jo Jo 


and l is defined by 
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We introduce the Hilbert space A = if*(0,1) x i7g(0, 1) x L 2 (0,1) equipped with the 
norm 

IMIa = ll^lx + U 3 II 2 + IHlI + \\ v 3x\\l + ||V5||1- 
It is clear that b is a bilinear and continuous form on A x A, and l is a linear and 
continuous form on A. Furthermore, there exists a positive constant Co such that 

b(v,v) = k\\v lx + v 3 \\l + Pi\\vi\\l + b\\v 3x \\l + (S 2 ((3 + t) + p 2 )\\v 3 \\l + p 3 \\v 5 \\l 

/* 1 PX p 1 PX 

+ 2 p 3 (/3 + t)6 / v 3 ( v 5 ds)dx + pl(/3 + t) / (/ v 5 ds) 2 dx 

Jo Jo Jo Jo 

> colMli- 

which implies that b is coercive. 

Therefore, using the Lax-Milgram theorem we conclude that the system (12.911 has a 
unique solution 

(vi,v 3 ,v 5 ) G (Hi(0,1) x H ( \(0, 1) x A 2 (0,1)), 
and we deduce from (I2.4I) - (I2.6I) the existence of v 2 G 0,1), tq G FTq( 0, 1), and 
v 6 G L 2 (0,1)) C L 2 (0,1)). 

Now, it remains to show that 

v\ G H 2 ( 0,1) n H}( 0,1), v 3 G H 2 ( 0,1) n H%( 0,1), v 5 G Hl( 0,1) and v 6 G -Hq(0 , l )- 
From (12.91) . we have 

~kv lxx = kv 3x - pivi + hi G L 2 ( 0,1). 

Consequently, it follows that 

vi e H 2 (o. i)nHl(o,i). 

Moreover, (|2.10p is also true for any ipi G C 1 ([0,1]). Hence, we have 


k / vi x tp\ X dx + k v 3 ipi x dx + p\ I vppidx = / hiipidx, 

Jo Jo Jo Jo 

for any </q € C 1 ([0,1]). Thus, using integration by parts we obtain 

vi x {l)ipi(l) - v\ x (0)(/?i(0) = 0, for all G C 1 ([0,1]). 
Therefore, v\ x {\) = vi x (0) = 0, and we deduce that 

vi G H 2 {0,1) r\Hl(0,l). 

Now, we substitute (|2.3p « into (12.311 .1 . we get 

bv 3xx = kv lx + kv 3 + 5twq - 5{fJ + t)v 6 + p 2 v 3 - h 2 G L 2 (0,1). 
Consequently, it follows that 

« 3 G^ 2 (0,1) 11^(0, !)• 

On the other hand, we get from (I2.3l) fi. 

V5x = tw 6 - (/3 + r)v e G L 2 (0,1), 


and we deduce that 

« 5 Gff 1 ( 0 ,l)n £ 2 ( 0 , 1 ). 

Similarly, from (12.31) we have 

V6 X = p 3 w 5 +5w 3x - dv 3x - p 3 v 5 <E L 2 (0,1)) which implies vq G Hq(0, 1), 
as u 6 (0) = u 6 (l) = 0. 

Finally, the operator I + A is surjective. ■ 
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Using Lemmas 12.11 and 12.21 we conclude that the operator A + B is the infinitesimal 
generator of a non-linear contraction Co-semi-group on the Hilbert space #. 

Finally, by applying the semi-group theory to m (see ID, we easily get the 
following well-posedness result. 

Theorem 2 . 1 . Assume that (Hi) and (A 2 ) are satisfied, then for all initial data 

(<po,fpi,^o,^i,do,qo) e (#*(0,1) n#*(o,i)) x #|(o,i) x (# 2 (o,i) n#o(o, 1)) 

x #o(0,1) x #*(0,1) x #o(0,1), 
the system il.8\) has a unique solution (ip, if, 6, q ) that verifies 

(ip, VO € c°(R+, (#2(0, 1 ) n #1(0, 1 )) X (#2(0, 1 ) n #^( 0 , 1 ))) 

n c^iL)-, #*(o, 1 ) x # 0 ( 0 , 1 )) n c 2 (m+,l 2 (o,i) xl 2 (o,i)), 

and 

(0,q)€ C°(R+,#i(0,l) x # 0 1 (0,1)) n C' 1 (M+,L 2 (0,1) xL 2 (0,1)). 

3. Stability results 

In this section, we state and prove a stability result for the nonlinear Timoshenko 
system (USD- For this purpose, we consider the following hypotheses: 

(Hi) : a : R + -A R_|_ is a differentiable and decreasing function. 

(H 2 )* : h : R + -A R+ is a continuous non-decreasing function such that h( 0) = 0 and 
there exists a continuous strictly increasing odd function ho € C([0, +oo)), continuously 
differentiable in a neighborhood of 0 and satisfying ho(0) = 0 

f h 0 (\(s)\) < |h(s)| < h^(\(s)\), for all \s\ < £, 

\ Ci|s| < \h(s)\ < C 2 1s|, for all |s| > e. 

where c* > 0 for i = 1 , 2 . 

Moreover, we define a function # by 

(3.1) H(x ) = y/xho(y/x) 

Thanks to Assumption (H 2 )*, # is of class C 1 and is strictly convex on ( 0 ,r 2 ], where 
r > 0 is a sufficiently small number. 

Remark 1. 


• We denote by c positive generic constant throughout this paper. 

• The hypothesis A\ implies that a(t) < c. 


We recall here the stability number defined by : 

/ Pi s/P2 Pis 

( W'T-fc 0 


t5 2 Pi 

bkp 3 


3.1. The case p = 0. 

In this part, we state and prove the decay results which are not necessarily of exponential 
or polynomial types. For this purpose, we establish several lemmas. We recall that the 
energy associated with the system (11.81) is defined by 

1 f 1 

(3.2) E(t) := - j ^ (pnp 2 t + p 2 fi 2 + bfil + k(ip x + if) 2 + p 3 6 2 + rq 2 ) dx. 
Throughout the rest of this paper we assume that conditions (Hi) and (H 2 )* hold. 
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Lemma 3.1. Let (ip , ip , 9 , q) be a solution of the system (11.81) . Then, the functional E 
satisfies 

(3.3) E'{t) = —/? f q 2 dx — aft) f tfth(ift)dx < 0. 

Jo Jo 

Proof. By multiplying the first fourth equations in (11.81) . respectively, by tp t , ip t , 6 and 
q, using the integration by parts with respect to x over ( 0 , 1 ), the boundary conditions 
(USDs and the hypotheses (Ai) and (A 2 )* , we obtain (13.31) . ■ 

Lemma 3.2. Let (ip, if>, 9, q) be a solution of the system (jl.8p . Then, the functional 

(3.4) R\(t) := - [ (pupipt + P2ipipt)dx, 

J 0 

verifies the following estimate 

(3.5) K[(t) < — Pi f <Pt’ dx — P 2 f iftdx + c f if 2 dx + k f {p x + f)) 2 dx 

Jo Jo Jo Jo 

+ | [ 6 2 dx + \ [ h 2 fif t )dx. 

z Jo 1 Jo 

Proof. By differentiating (13.41) and using the first and second equations of (11.81) , we get 

K[(t) = ~ Pi Ptdx-p 2 / V’t dx — / k(<p x + i/)) x <pdx — / (bif xx - k(<p x + ip) 

Jo Jo Jo Jo 

— 59 x — a{t)hff>t))f>dx. 

Integrating by parts and using the boundary conditions (|1.8jh . we have 

K[ (t) = —p\ f tp 2 dx — P 2 f f> 2 dx + b f 'f 2 dx + f k(<p x + if) 2 dx 

Jo Jo Jo Jo 

—6 / 0f> x dx + / a(t)h('il’t)'f>dx. 

Jo Jo 

Applying Young’s inequality, we obtain (13.5D . ■ 


Lemma 3.3. Let (p,if,0,q) be a solution of the system (11.81) . Then, the functional 

(3.6) K 2 (t) \= p 2 / ifnftdx — p 2 / <p t wdx - 5 t / fjqdx, 

Jo Jo Jo 

satisfies, for any e > 0 

(3.7) K' 2 (t) < — (b — 2ce) f i/j 2 dx + c(f if 2 dx+ f q 2 dx+ f h 2 (ipt)dx) 

Jo Jo Jo Jo 


+ Pis' / <Ptdx, 

Jo 

where w is the solution of the problem 

^XX — f’xi 

u>(0) = u;(l) = 0. 


(3.8) 
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Proof. By differentiation of (13.61) and the use of the first, second and fourth equations 
of dll), we get 

p 1 pi pi pi pi 

K' 2 (t) = P 2 fi’ 2 dx+b / if xx ifdx—k / (tfix+ip^dx—d / 9 x ifdx—a{t) / iph(ipt)dx 

Jo Jo Jo Jo Jo 

pi pi pi pi pi 

+k / (<£> 3 , + fi>) x wdx + pi tptwtdx — t 5 iftqdx + 5/3 ipqdx + 5 6 x ifdx. 

Jo Jo Jo Jo Jo 

Integrating by parts the last equality, using (13.81) and the boundary conditions (II. 81b . 
we have 

pi pi pi pi pi 

K' 2 (t) = p 2 / if 2 dx — b / '(/j x dx — k / fi> 2 dx + k / w 2 dx — aft) / iph(ipt)dx 

Jo Jo Jo Jo Jo 

+Pi / PtWtdx — t 5 iftqdx + 5(3 tfqdx. 

Jo Jo Jo 

By a simple calculation, we easily deduce that the function w satisfies the following 


estimates 




f 1 

/*! 

(3.9) 

/ w 2 dx < 

/ 'ffidx, 


Jo J 

'o 


f 1 

ri ^ 

(3.10) 

/ wJdx < c 

Jo 

Jo ^ dX 


Thanks to Young’s and Poincare’s inequalities and (13.9D - (13. 1Q[) . we conclude that 

(3.11) K 2 (t) < p 2 ['ipfdx — b ffi 2 dx + y- f w 2 dx + p\e [ <p 2 dx 

Jo Jo 4e Jo Jo 

+ t5s / i(> 2 dx + 


-. 2 ^ , r& r 1 
Jo 


4e 


q 2 dx + c p £ I ip x dx + J q 2 dx 


f 1 c 2 f 1 

+ £c p / tpldx + — / h 2 (ift)dx. 

Jo Jo 

Therefore, we obtain (13.71) . 


Lemma 3.4. Let (ip,fi,6,q) be a solution of the system (|1.8D . Then, the functional 

1 px 

q{ 6(t,y)dy)dx , 

■Jo 

satisfies 


(3.12) 


K 3 (t) := -rp 3 


(3.13) K' 3 (t)<~y- f 9 2 dx + c( f q 2 dx+ f fi> 2 dx 

2 Jo \J o Jo 

Proof. By differentiation of (13.121) and the use of the third and fourth equations of (11.81) , 
we get 



K’ 3 (t) 



0(t,y)dy)dx + p 3 


q x (t,y)dy)dx + t5 


f o x ( r 

Jo Jo 

pi px 

/ q{ 

Jo Jo 


0{t,y)dy)dx 

iftx{t,y)dy)dx. 
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By integrating the above equality over (0,1) and using the boundary conditions (II.8L 
(note also that j 0 * ddx = 0), we have 

/*1 px rl /*1 /»1 

K' 3 (t) = ps(3 / q( 9(t,y)dy)dx - p 3 / 6 2 dx + r / q 2 dx + t5 / qip t dx. 

Jo Jo Jo Jo Jo 

Applying again Young’s inequality and the fact that 

rl rx rl 

/ (/ d(t,y)dy) 2 dx < c / 0 2 dx , 

Jo Jo Jo 

we arrive at (J3.13I) . ■ 


Lemma 3.5. Let 0,q) be a solution of the system (II. 8jl . Then, the functional 

brpi rl 


(3.14) K 4 (t) := ^ J iftiVx + ip)dx + J iptijxdx 


U2 

K ■! 0 
)T 

5k y b k jo 


-jjrijf ~ y)/ e vtdx + q{}Px + ^)dx, 


satisfies 


(3.15) ^(t) < —(r — 2e!) f(ip x + i/j) 2 dx + C f I 'i/j 2 dx+ f q 2 dx+ fh 2 ( K 'ip t )da 

Jo \J 0 4o do 

/ 0x(<Px + ip)dx, 

Jo 


bp3 

Spi 


( pl )( p2 Pl ) T62pl 

kp 3 b k bkp 3 


with C = 2max(f + ±, (^(f - f )) 2 (J^ + £), and e x > 0. 

Proof. By differentiation of (13.141) . using (11.81) and integration over (0,1), we get 

K' A (t) = ^ [ {bf’ xx - k(<p x + if)~ 56 x - a(t)h(pp t )){.Tx + ^)dx 

4 Jo 

[ 1pt(Vx+Tp)t ’ 

Jo 


TP 2 

k 


) t dx + ^2 y (Tx+ i>)xVx + Vti’txdx 

7 /* 1 

- Tt(t--T-) (-(qx + ty x t)Tt + 0{Vx + fj)x)dx 

dk 0 k Jo 

+ Jk^T ~ J 0 + 9x ^ ipx + ^) + + J>)t)dx. 

By integration over (0,1) and using the boundary conditions (II.8L . we have 


-^ 4 ( 0 = ~T ( tp x + i/j) 2 dx + 




5k b k jo 

^(y ~ J 0 q ^ x + ^ dx ~~p J 0 ot{f)h{i> t )^ x + if)dx 


+ 


V 3 

<5pi 


/ Pi wP2 Pi^ rd 2 pi 

kp 3 b k bkp 3 j jo 


^(Px + f))dx. 


Applying Young’s inequality, we obtain (13.151) . 


Next, we define a Lyapunov functional AT and show that it is equivalent to the energy 
functional E. 
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Lemma 3.6. Let (ip,if, 6 , q) be a solution of the system (11.81) . Then, the functional 

(3.16) K(t) := NE(t) + K x + N 2 K 2 + N 3 K 3 + N 4 K 4 , 

where N is sufficiently large, N 4 and N 2 are positive real numbers to be chosen properly, 
satisfies 

(3.17) c\E(t) < K(t ) < c 2 E(t), 
for ci and c 2 two positive constants and 

(3.18) K'(t) < — (pi — N 2 p 4 £) f pfdx — p 2 f iffdx — (N 2 (b — 2ce) — c) f iffdx 

Jo Jo Jo 

- [ (AT 4 (t- 2s 4 ) - k))(<p x + iffdx- (—^ - [ 9 2 dx 

Jo 1 1 Jo 


— (N(3 — cN 2 — cN 3 — cN 4 ) f q 2 dx + c f (iff + h 2 (if t ))dx 

Jo Jo 


+ N * b ir 

dpi 


( _ _ Pi \ _ t ^ 2 pi 

kp 3 b k bkp 3 


-1 


O x (<Px + if)dx. 


Proof. From Lemmas 13.21 to [331 we find 

\K(t) - NE(t)\ < pi f \<pip t \dx + (p 2 + N 2 ) [ \ifif t \dx + N 2 p 4 f \<p t w\da 

Jo Jo Jo 

rl rl rx 

+ N 2 t5 / \ifq\dx + rp 3 \q( 6(t,y)dy)\dx. 

Jo Jo Jo 

Applying Young, Poincare and Cauchy-Schwartz inequalities and the fact that 
Tx — ^(Px + if f + 2 if 2 < 2(<p x + if) 2 + 2 ciff., 
we obtain (13.171) . and therefore we get 

K(t) ~ E(t). 


For to prove (13.181) , it suffices to differentiate (13.161) and use lemmas 13.1113.51 This ends 
the proof of the lemma. ■ 


Theorem 3.1. Let us suppose that 

P = 


, Pi wP2 Pi x t5 2 pi 


kp 3 b k bkp 3 


= 0 . 


Then there exist positive constants k\, k 2 , k 3 and £q such that the energy E(t) associated 
with (USD satisfies 

(3.19) E(t) < k 3 Hf l J a(s ) ds + k^ , for all t > 0, 


where 


ri i t 

Hi(t) = / -—-ds, H 2 (t) =tH' (e 0 t). 

Jt H 2 (s) 


Here H 4 is a strictly decreasing and convex function on (0,1], with limfl’i(i) = +oo. 
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Proof. The estimate (13.181) . with p = 0, takes the form 

K'(t)< - {pi-N 2 pi£) [ (fitdx - p 2 [ iff dx — (JV 2 (6 — 2ce) — c) [ if; 

Jo Jo Jo 

- f (N 4 (t - 2ei) - fc))(<Az + if) 2 dx - (—^ - |) / d 2 dx 
Jo z z Jo 


dx 


- (N/3 - cN 2 - cN 3 - cN 4 ) f q 2 dx + c f (i/j 2 + K 2 {ift))dx. 

Jo Jo 

Now, we choose the constants in the above estimate as follows: first e and e i are such 
that 


1 


£ = 


2 N 2 


and £i < —. 

2 


After that, we choose N, N 2 , N 3 and N 4 sufficiently large such that N 2 > -ff, N 3 > —, 
N 4 > r _ k 2 £ 1 and N > ^ + T _ 2 El )• Then, we deduce that 

(3.20) K\t) < —dE(t) + c f {■$ + h 2 (if t ))dx, 

Jo 


where d = min(pi — N 2 p \£, p 2 ,N 2 (b — 2ce) — c, N 4 (t — 2e\) — k, N ?f 3 — | ,Nj3 — cN 2 — 
cN 3 — cN 4 ). 

First case: Let ho be a linear function over [0, e]. The hypothesis (A 2 )* implies that 

c' x |s| < \h(s)\ < c 2 1s|, for all s € R. 

Consequently, by multiplying inequality (13.201) by aft), we obtain 

(3.21) a(t)K'ft) < —daft)Eft) + c aft) f (ip 2 + h 2 (ip t ))dx, 

Jo 

r 1 1 

< —daft)Eft) + caft) / (— \ip t h(ip t )\ + ^l^tK^Ddx, 

Jo c x 

< — da ft) Eft) + coaft) / ipth(ipt)dx = —da(t)E(t) — cqE 1 ft), 

Jo 

where Co = c(4- + c 2 ). 

C 1 

Using now hypotesis (Ai), this yields 

(3.22) (aK + coE)'ft) < a(t)K'(t) + c 0 E'(t) < -da(t)E(t). 

We integrate the inequality (13.221) and use the fact that aK + cqE ~ E, we 
obtain for some &, c > 0, 


(3.23) 


Eft) < fcex 



a(s)cZs). 


Finally, by a simple computation we get (13.191) . 

Second case: Let ho be a non-linear function over [0, e]. We assume that max(r,ho(r)) < e, 
where r is defined in the hypothesis (A 2 )*. 

Let £\ = min(r,ho(r)), we deduce from the hypothesis (A 2 )* that 


Mfilu < LP| S | < | ftW | < < IeMl, 

£ |s| |s| £1 
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for all s satisfying £1 < |s| < e. 

Then, the estimates in hypothesis ( Ai)* become 


(3.24) 


f ho(\s\) < |h(s)| < h 0 \\s\), for all |s| < £ 1 , 

\ c x |s| < \h(s)\ < c 2 |s|, for all |s| > £ 1 , 


and we have 


(3.25) s 2 + h 2 (s) < 2H~ 1 (sh(s)). 

To estimate the last term of (13.201) . we consider the following partition of (0.1): 

^1 = {x € (0,1); IV’tl < £ 1 }, 0 2 = {x € (0,1); \i/j t \ > £ 1 }. 

Then, we obtain 

(3.26) ipthfyt) < H{r 2 ) and < x 2 on Q]. 


Now, we apply Jensen’s inequality to the following term 
I(t) := |^| ip t h{^t)dx, 


and we infer that 

(3.27) H~\l(t)) > c [ H~\ip t h(ilH))dx. 

Using (13.24j) . (I3.25|) and (|3.27|l . then the right-hand side of (I3.20P multiplied by 
a(t) becomes 

a(t) f (i/j 2 + h 2 (i/j t ))dx = a(t ) f (i/j 2 + h 2 (i/j t ))dx + a(t) f (if ) 2 + /i 2 (V>t))dx, 

Jo J fii J n 2 

< 2 a{t) f H~ l (^th(ilH))dx 

J f2i 


“h Ot 


( t ) f {\^t\^\K^ t )\+ ^t\\h{^ t )\)dx, 
J ^2 ^1 


< ca(t)H 1 (I (t)) + a(t)c iptH^dx, 

Jo 

< ca(t)H~ 1 (I(t)) — cE'(t). 
Consequently, the estimate (13.201) gives 

(3.28) i?o(t) < —da{t)E(t) + ca{t)H" 1 {I(t)) , 


where Ro = aK + cE. 

On the one hand, for £q < r 2 , using (13.281) . H' > 0 and H >0 over (0, r 2 ] and 
E' < 0 the functional R\ defined by 

Ri(t) := H'^£ 0 W^jR 0 (t)+c 0 E(t), 

is equivalent to E(t). 

On the other hand, using the fact that g o H" ( g o ~Ej§j)Ro(t) < 0 and (13.281) . 
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we conclude that 

(3.29) R[(t) = e 0 ^-H"(E 0 ^)Ro(t) + H\e 0 ^)R' 0 (t) + coE\t) 

< -da(t)E(t)H'(e 0 |||) + ca{t)H\e 0 ^)H-\l(t)) + c 0 E'(t). 

Our goal now is to estimate the second term in the right-hand side of (13.291) . 
For that purpose, we introduce the convex conjugate H* of H defined by 

(3.30) H*(s) = s(H’)-\s) - F((FT , ) _1 (s) for s € (0, H'(r 2 )), 
and H* satisfies the following Young inequality: 

(3.31) AB < H*(A) + H(B) for A € (0, H (r 2 )), B € (0,r 2 ). 

Now, taking A = H'(e o-g^y) and B = t)), we obtain 

K(t) < + 

+ca(t)H + c 0 E\t) 

< -daWW'ie 0 ||) + ceo || a(i) ^( £o ||, 

-ca{t)H(e 0 E^~) + ca{t)I(t) + c 0 E r (t) 

E(U) 

< -ia(t)E{t)H'(e :„§j|) + c£o|||) - cE'(t) + n E'(t). 

With a suitable choice of £o and co, we deduce from the last inequality that 

(3.32) R[(t) < ~(dE( 0) - aro)a(t)||)/f'(£(,|||) < -to(t)H 2 (|||), 
where k = dE{ 0) — ceq > 0 and ^(s) = sH' (eqs). 


Since E{t) ~ i?i(t), then there exist ai and 02 such that 
ai-Ri(t) < E(t) < d2Ri(t). 

We set now R(t) = ■ It is clear that R(t) ~ E(t). We use the fact that 

H 2 (t), H 2 (t) > 0 over (0,1] (this is due to the fact that H is strictly convex on 
(0,r 2 ]) and we deduce from (|3.44l) that 

R'(t) < —kia(t)H 2 (R(t)), for all t € M+, 

with k\ > 0 . 

By integrating the last inequality, we obtain 


Hi(R(t)) > H!(R(0)) + h f a(s) ds. 

Jo 

Finally, using the fact that H^ 1 is decreasing (because H\ is also), we have 

R(t) < H- { J* a(s) ds + &) 2 ^ , with k 2 > 0 . 

Taking into account that E(t) ~ R(t), we deduce (|3.19l) . 
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3.1.1. Examples. In the following, we will apply the inequality (|3.19l) on some examples 
in order to show explicit stability results in term of asymptotic profiles in time. For 
that, we choose the function H strictly convex near zero. 

Example 1. 

Let h be a function that satisfies 

C 3 min(|s|, |s| p ) < \h(s)\ < C 4 max(|s|, |s| p ), 
with some C3, C4 > 0 and p > 1 . 

For /iq(s) = cs p , hypothesis (A 2 )* is verihed. Then H(s) = cs^~. 


Therefore, we distinguish the following two cases: 

• If p=l, we have ho is linear and H 2 (s) = cs, H\(s) = — and FI 1 _ 1 (t) = 
exp(— ct). 

Applying (13.191) of Theorem 13.11 we conclude that 

E(t) < k^ exp(— c(k\ / a{s)ds + k 2 )). 


and 


p — 1 p_^ 

lip > 1 ; this implies that ho is nonlinear and we have H 2 (s) = c^^-e 0 2 s~ 2 ~ 


Therefore, 


2 p— 1 


HM = I, = loh) ~ WT) 


t 2 , with 5 = c 


p+ 1 




Example 2. 


Using again (13.191) . we obtain 

E(t) < H^ 1 (k\ f a(s) ds + k 2 ) = (S^——(hi [ a(s) ds + fe) + 1)~ 

Jo 2 Jo 

Let ho(s) = exp(—i), this yields H(s) = • v /sexp(— j=) and 


2 

p-1 


H 2 {s) = + -h) eX p(--L=). 

2y/s 0 2e 0 v £ o s 


Moreover, we have 

Hi(t) = 

< 


I! ( 3Tj) Hp( ii |ds 


< c 


1 V- 3 

\ 2-v/so 

+ 

2 \/en 


r ex P 
Vs 

1 1 



f 1 11 1 1 

/ 0 - exp(— =)ds = cexp(^=) - cexp(^). 

Jt 2s^/So s \/ £ 0 s V e ot v^o 


Then, 


t < £ 0 l (in 


H\{t) + cexp(^=j)\ \ 
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h; 


Replacing t by H 1 1 ^k\ J a(s) ds + k 2 ^ in the last inequality, we find 

( ( 1 \ \ ~ 2 

' / a(s) ds + k 2 + cexp(——) 

Jo V^o 


rr ‘H a(s) ds + k 2 ) < e. 
Therefore, 

-E(t) < ^3^0 1 

Example 3. 


-1 

0 


/ / y* 1 \ \ 

A'i / a(s) ds + ^2 + cexp(——) 

Jo v^o 


-2 


In 


V 


V 


/ 


Let do(s) = 7 exp(—4j). Following the same steps in exemple 2 we find that the 
energy of (j 1 . 8 f) satisfies 

/ ki ft a(s) ds + k 2 + cexp(i) \ 


E(t) < e ln(- 


Example 4. 


Let Lo('S) = | exp(—|(ln s) 2 ). Then, we have H(s) = exp(—|(ln s) 2 ), 


H 2 (s) = - 


s ^ \ 4 

_ _1 In eps 


2 eo 

4 4 s 


ex P(-l( ln£ os) 2 ) and idi(t) = fl -2 I ^exp(i(lne 0 s) 2 ). 


As lim —— , u = 0, then the function s i-A ,, is bounded on (0, ll, and 

s — >0 (lnOos )) 2 v ’ J ’ 

we infer that 


Hi(t) < c] - 
Hence, we have 


1 ln£o£ eX p(l(i n s ) 2 ) ds = eX p(^(lne 0 i) 2 ) - exp( j(lne 0 ) 2 ) • 
t 1 £o« 4 4 ._ 4 _, 


t < — exp ( -2 (ln(idi(f)) + ci )2 
£0 


Replacing t by H 


1 1 (ki J a{s) 


ds + k 2 I in the last inequality, we find 


Eft) < ksH 1 1 ^k\ J a(s) ds + k^j = — exp ^—2 ^ln k\ J a(s) ds + k 2 + J . 
3.2. The case p / 0 and aft) = 1. 

This section is devoted to the statement and the proof of the stability result for the 
system (j 1.8 1 ) when p / 0 and aft) = 1. 

We have the following theorem. 


Theorem 3.2. Let us suppose that conditions (Hi) and (A 2 )* hold, then for 

/ 0 , 


, P\ U P2 Pi^ t5 2 pi 

( r ~ r— )(-r - 


kp3 b k bkp^ 


P = 
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the energy solution of (11.811 satisfies 


(3.33) E(t)<Hfi\-), 


where 


H 2 (t) = tH (eo t) with lim H 2 (t) = 0. 


t—>0 


Proof. Let (ip,ip,9,q) be a solution of the system (11.81) . First, we define 
1 f 1 

(3.34) Eft) ■= - {piPt + P^t + H’l + KPx + ip ) 2 + P 3 O 2 + rq 2 ) dx, 

1 Jo 

and 

If 1 

(3.35) Eft) := - jI {pi<p 2 t + p 2 ip 2 t + biptx + Hvtx + 'fitf + P3&t + T Qt) dx. 

>atisfies 

'(t) = -fi [ q 2 dx — [ ip t h(ip t )di 

Jo Jo 


Then, the functional E satisfies 
(3.36) E'(t) = -fi f q 2 dx — 


ix < 0. 


Analogously, the functional E satisfies 

(3.37) E’(t) = -/3 f qfdx - [ 'ip 2 t ti (fifid, 

Jo Jo 


x < 0. 


Using the results in Subsection 13.II (recall the expressions of the functionals K\, ...,K 4 ) 
we have the following Lemma. 

Lemma 3.7. Let (ip,ip,6,q) be a solution of the system (|1.8D . Then, the functional 

(3.38) L(t) ■= N(E(t) + E(t)) + K x + N 2 K 2 + N 3 K 3 + N 4 K 4 , 

satisfies 


(3.39) 


L'(t) < 


d'(t ) + c f (ip 2 + h 2 (ip t ))dx , 
Jo 


for N large enough and d >0. 

Proof. By differentiation of (I3.38p . and using (13.181) and Young’s inequality, we obtain 


(3.40) L'(t ) < —dE(t) + c f (ip 2 + h 2 (ip t ))dx + c f (9 2 + (<p x + ip) 2 )dx 

Jo Jo 

r l rl 

- 2 ' 7 ™ AT / ip 2 t h'(ipt)d. 

Jo 


—N/3 / q 2 dx — N 


x. 


10 


Now, from (ll.SD .i . we deduce that 


J 0 2 dx < c( [ q 2 dx + [ q 2 dx 


1 0 


>0 


Consequently, we get 

(3.41) L'(t) < —d!E(t) + c f (ip 2 + h 2 (ipt))dx — (fiN — c) f q 2 dx 

Jo Jo 

~ N [ iplth'(fi>t)dx. 

Jo 






22 


AYADI, BCHATNIA, HAMOUDA, AND MESSAOUDI 


where d! = d — c > 0 and d is the same constant that appears in (13.201) . Finally, we 
choose N large enough and using the monotonie of the function h we arrive at (|3.39|) . 


Now, using the following partion of (0,1) defined in Subsection 13.11 the right-hand 
side of (13.391) becomes 

[ + h 2 (i/j t ))dx = f {i$ + h 2 (il> t ))dx + [ (V| + h 2 (i> t ))dx. 

Jo Jn i Jfio 


/o Jn i 

Now, the estimates (|3.24D - (I3.27I) imply that 


Jo 


+ h 2 (i/j t ))dx <2 H 1 ('ip t h(ip t ))dx 


’ f2i 


+ 

< cH 


[ m\m)\+ 

Jn 2 ci 

1 (I(t)) + c[ il) t h(il)t)dx. 

Jo 


Consequently, 


L'(t) < —d'E(t) + cH 1 (I(t)) + c f i/jth(i/jt)dx + c/5 f q 2 dx 

Jo Jo 

< -d'E(t ) + cH~ l (I(t )) - cE'{t). 

Hence, we deduce that 

(3.42) (L + cE)'{t) < -d'E{t) + cH~ l {I{t)). 

We then define 

Ri(t) ■= H'(e 0 + cE)(t) + co E(t), 

which verifies 

(3.43) R[(t) < -d 1 E{t)H\e G ^) + cH\e 0 ^)H-\l(t)) + eE\tl 
as we have e 0 ^yH" (e 0 §^)Ro(t) < 0. 

We recall the definition of the convex conjugate H* of H, given by (I3.30p . which satisfies 
the following Young inequality: 

AB < H*(A) + H(B) for A € (0,^'(r 2 )), B € (0,r 2 ). 

With the same choice of A and B as in (13.311) . we obtain 

With a suitable choice of Co and e, we deduce from the above inequality that 

(3.44) R[(t) < ~(dE( 0) - “o)|^ff'(Eo|^) < 
where k = dE{ 0) — ceo > 0 and i?2(s) = sH (eq(s)). 
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Finally, we have 


R[(t) < -hH2 (||y), for all t € M+, 


with fci > 0, which yields 

< J‘ < _(fl l(t ) - fil (0)) < *(0). 

Then, we easily deduce that 


£(t) ^ i?i(0) 

ffa( l(0) ) 5 


Thus, 

This concludes the proof of Theorem 13.21 




3.2.1. Examples. 

p -\-1 

Example 1: Let /io(s) = cs p . Then H(s) = cs~ 2 ~. 

Therefore, we distinguish the following two cases: 

• If p=l, we have ho is linear and H^it) = cs. 
Applying (13.331) of Theorem 13.21 we conclude that 

m < f 


P-1 


• If p > 1; this implies that ho is nonlinear and we have H 2 (s) = cs 2 . Therefore, 
Using (13.331) . we obtain 

E(t) < ct~p 31 . 

Examples 2: Let h be given by h(x) = exp(—-^) and we choose ho(x) = 1+ jf exp(—^r), 

gxp ( _-_^ 

we obtain iL(x) = ±±^exp(-±) and iL 2 (x) = — s sJ 0X ■ 

Then, we use the following property : 


and we deduce that 


lim exp(- )H 2 (x) = + 00 , 

£->■0+ £qX 


exp(-3-) < H 2 (x). 

£ 0 X 


We infer that there exists xq > 0 such that, 

exp(-—) < H 2 (x ) on (0,x 0 ]. 

EqX 

Consequently, the energy of the solution of (11.8j) satisfies the estimate 


E(t) < c(ln(i)) 1 . 
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